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We use molecular dynamics simulations to investigate the microscopic and macroscopic response 
of model polymer networks to uniaxial elongations. By studying networks with strands lengths 
ranging from N 3 = 20 to 200 we cover the full crossover from cross-link to entanglement dominated 
behavior. Our results support a recent version of the tube model which accounts for the different 
strain dependence of chain localization due to chemical cross-links and entanglements. 
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Cross-linking a melt of linear precursor chains leads 
to a polymer network which macroscopically behaves as 
a (viscoelastic) solid and which is microscopically char- 
acterized by a complex, quenched, random connectivity 
and topology Q. Over the past sixty years a large va- 
riety of theories of rubber elasticity has been put for- 
ward HHQt, and there is a corresponding body of rhe- 
ological literature devoted to comparing and testing the 
proposed stress-strain relations, see e.g. jE 0|. Only 
recently, neutron scatte ring exp eriments and com- 

puter simulations HE1 llilliiLUI have begun to provide 
detailed microscopic information. Attempts to quantita- 
tively correlate the microscopic and macroscopic response 
to strain were so far restricted to idealized model polymer 
networks with diamond lattice connectivity [llj, LL2J ■ 

In this Letter we report a comprehensive set of results 
from computer simulations of randomly end- linked model 
polymer networks 0, 0] under elongational strain. In 
our data analysis we follow the logic of most statisti- 
cal mechanical theories of rubber elasticity 0, 0| i-e. 
we relate the localization of different parts of a poly- 
mer network to strain-induced microscopic deformations 
and the macroscopic elastic resp onse 0, Q] . We focus 
on theories EUl El El El El El based on Edwards' 
tube model |20j, because the underlying ideas are con- 
ceptually relatively simple, (almost) completely worked 
out, and closely related to most modern theories of poly- 
mer rheology 00]. Nevertheless, a key ingredient of 
the model, the strain and stand length dependence of 
the phenomenological tube diameter, remains controver- 
sial [H Our data support a recent generalization 
of the tube model Ell which accounts for the different 
character of chain localization by cross- links [lfil and 
entanglements |2.llsUlfl|. 

We used extensive Molecular Dynamics (MD) simula- 
tions to study the behavior of bead-spring polymer melts 
under uni-axial, volume-conserving elongation. The 
polymer model has two interactions that represent bonds 
and excluded volume, respectively. Interaction parame- 
ters are chosen to ensure conservation of the topological 



state. The details of the simulation methodology can be 
found in the literature, see e.g. El El- The present 
Letter is based on simulations of end-linked networks of 
M x N s = 5000 x 20, 1000 x 35, 2500 x 100, and 3000 x 200 
where M denotes the number of strands and N s strand 
length. The state of the networks can be accurately char- 
acterized: the gel faction is > 99%, > 91% of the network 
is elastically active, and the fraction of four-functional 
crosslinkers is > 78%, except for 62% in the case of the 
3000 x 200. Each network was simulated at several elon- 
gations (Ax =A|| = A, X y =X Z =Aj_ =l/VX). The net- 
works were successively strained. After each strain incre- 
ment the network was equilibrated, and configurations 
were samped for more than 30 (2) Rouse times of the 
strands for the two short (long) strand networks, respec- 
tively. To reduce finite chain length effects the maximal 
elongations were limited to A = 2 (A = 4) for short (long) 
strand networks, respectively. Elastic properties were ob- 
tained by sampling the deviatoric part of the microscopic 

virial tensor defined as o a p — Fij,a r ij,0^ /V ■> where 

the sum is over all pairs i,j of interacting beads, a, f3 are 
Cartesian indices, and F, r and V denote forces, separa- 
tions and the volume, respectively. The normal tension 
is defined as <tt = &xx ~ {pyy + °2z)/2. In the following 
we present our simulation results together with a brief 
outline of the theoretical background. 

With the tube model Edwards proposed a 
tractable simplification of the complicated many-body 
problem of a randomly cross-linked and entangled poly- 
mer network. The key idea is that vulcanization leads 
to a permanent localization of the precursor chains (or, 
equivalently, long, randomly chosen paths through the 
network) in tube-like regions along the coarse-grained 
chain contours. In the following d(X a ) denotes the cc'th 
Cartesian component of the tube diameter, which de- 
pends only on the corresponding strain component A a , 
hence each elongation provides two tube diameters rf(A||) 
and d(X±). The zero-strain tube diameter d(X a = 1) is 
abbreviated d. The simplest measure of the tube diam- 
eter is the width, dx{X a ), of the thermal fluctuations 
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FIG. 1: (color online) Strain-dependent localization of cross- 
links for N a = 20 (magenta A), N s = 35 (red o), N s = 100 
(green □), and N a = 200 (blue o). (a) Cross-link mean- 
square displacements as a function of time for t he un strained 
networks. Lines are fits of 2d^ r (A Q )[l — exp(— tfrx)] which 
we use to extract the tube diameters, (b) Rescaled cross-link 
mean-square displacements for all {N a , \ a } in comparison to 
the fit function, (c) Strain dependence of the tube diameter. 
The lines are fits of Eq. (d) Strand length dependence 

of crosslinker confinement dx (blue *), dx,A component (red 
+), and dx.B component (green x). 




FIG. 2: (color online) Length scale dependent microscopic de- 
formation response (symbols as in Fig. |J). (a) Scaling plot of 
the degree of affmeness of the microscopic deformations using 
the tube diameters from Fig. The lines are /a (dashed) 
and /s (solid), (b) Parallel and perpendicular microscopic de- 
formations for the 2500 x 100 system for strain A = 1.5, 2, 3, 4 
(symbols). Lines are fit of Eq. 0. (c) Strain dependence of 
the tube diameters. Symbols indicate the error bars and the 
range of \ a covered by the data, (d) Strand length depen- 
dence of tube confinement dr (blue *), dr,A component (red 
+), and dr,B component (green x). 



of chemical cross-links around their average positions 
(Fig.^J. Our simulation runs are long enough to allow for 
a reliable extraction of dx from the components of the 
mean-square displacements gi, a {t) = ([r a (t] — r a (0)] ), 
where the average is restricted to four-functional cross- 
links (Figs.QJi and b). 

The key issue is the non-trivial dependence of the 
extracted tube diameters on strain and strand length 
(Fig-K and d). From a theoretical point of view, the sit- 
uation is relatively clear in the hypothetical case of "phan- 
tom networks" which consist of non-interacting Gaussian 
polymer chains [2 J| . The corresponding tube theory by 
Warner and Edwards (WE) uses an isotropic, strain 
independent tube diameter dx,A(^a) = dx,A, where 
dx,A ~ b^/W s is of the order of the root-mean-square 
extension of the network strands. The description of 
polymer networks as phantom networks becomes inap- 
propriate, if the length of the network strands approaches 
the melt entanglement length, N e fl4l l2l| . In partic- 
ular, the tube diameter should become independent of 
strand length in the limit of very long strands where 
dx,B ~ by/Ne. The strain dependence of entanglement 
dominated confinement is a subtle point. Empirical evi- 



dence and theoretical arguments sup- 

port the tube theories by Heinrich and Straube (HS) [2] 
and by Rubinstein and Panyukov (RP) ^3] which pre- 
dict an anisotropic, strain dependent tube diameter of 

dx,B(^a) = V\^dx,B- 

Instead of combining an entanglement tube model 
with a classical theory of rubber elasticity in an ad-hoc 
fashion, the recently introduced double tube theory by 
Mergell and Everaers (ME) jl^l accounts for the simul- 
taneous presence of both effects with the help two cor- 
related and additive confining potentials of the WE and 
the HS-RP type, respectively. The ME theory predicts 



4(A Q ) d\ Xld% ' v ; 

with an effective zero-strain tube diameter d^ 4 = d^ 4 + 
dg 4 . As indicated by the solid lines in Fig. QJ;, the in- 
dividual data sets are well described by the functional 
form of Eq. QJ. Moreover, the strand length dependence 
of the fitted cross-link and entanglement tube diameters 
dx,A and dx,B is in good agreement with the arguments 
underlying the double tube theory (Fig. QJl). It is worth 
noting the introduction of crosslinkers cause extra topo- 
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logical constraints, hence the observed narrowing of the 
entanglement tube for N s = 20. 

In the second part of our analysis, we consider mi- 
croscopic deformations. In rubber- like materials, macro- 
scopic deformations affect distances between neighboring 
monomers only weakly, while distances between distant 
monomers change affinely with the macroscopic strain. 
Following the logic of the tube model, we focus on the 
length scale dependent deformations of non-reversal ran- 
dom walk paths through the network. As a probe for 
microscopic deformations we utilize the mean-square dis- 
tances A a ) between pairs of beads as a function of 
their chemical distance n. In a Gaussian theory r 2 (n, A Q ) 
fully specifies the microscopic conformations. In the ab- 
sence of strain r^(n) = b 2 n, where b 2 is one Cartesian 
component of the mean-square segment length for a pre- 
cursor chain. End-linking is known not to ch ang e the 
segment length compared to the precursor melt . The 
crossover to afflne deformations can be characterized by 
a dimensionless function lial 



/ = 



X a ) - r 2 a {n) 



(2) 



with < / < 1. According to the tube model the 
crossover length should be of the order of dx, i.e. / 

„2 , 



should be a function of yx — r'a{n)/[2d x }. Fig. Eh ver- 
ifies the interdependence of localization and microscopic 
deformations, which is an essential element of the tube 
model. Note that the plot contains all available data sets 
({N s ,X a }), 

Theories based on the tube model make explicit pre- 
dictions for the functional form of / (see the lines in 
Fig. Hi): f A (y) = 1 + (exp(-y) - l)Jy Q and f B {y) = 
l + 0.5exp(— y) + 1.5(exp(— y)— l)/y [li| were obtained in 
the limits of cross-link and entanglement dominated con- 
finement, respectively. For the general case, the double 
tube model [l^J predicts 



f{y) = f A {y) + 



d 4 

a T,A 



A 4 \2 
U T,B a 



^T,A 



lf B (y) - f A (y)] , (3) 



where y = r^,(n)/[2(i^(A a )] with d T (X a ) given by Eq. JQ. 

We determined dx : A and g?t,b for each system by fit- 
ting Eq. (pj to the sampled r 2 (n, X a ) for all values of X a 
simultaneously |3fl |. A typical result is shown in Fig. Eh>- 
The effective strain-dependent tube diameters are shown 
in FigEh, the two components dx^A/ b{N s ) are plotted in 
FigEJi. As expected from the scaling plot Fig. Eh, there 
is qualitative agreement between the tube diameters dx 
extracted from the cross-link fluctuations and c?t inferred 
from the analysis of the microsopic deformations. In par- 
ticular, we observe a similar crossover to entanglement 
dominated confinement for long strands. Not surpris- 
ingly, there is no quantitative agreement between the two 
measures of the tube diameter. Tube models of rubber 
elasticity do not distinguish between middle monomers 



of network strands (which are free to slide along the en- 
tanglement tube) and cross-linkers (which are not) fl(i| . 
In fact, the theories discussed in this Letter employ a 
simple harmonic localization potential which suppresses 
reptation-like motion. 

In the third part of our analysis we consider the macro- 
scopic elastic response of our networks. Fig. El shows the 
strain-dependence of the sampled normal tensions in the 
form of a standard Mooney-Rivlin plot 0. This repre- 
sentation is commonly used to emphasize deviations from 
the classical stress-strain relation or (A) oc (A 2 — A -1 ). 
Similarly to experiments p], we find that these devia- 
tions are more pronounced for entanglement dominated 
systems and that the shear modulus decreases with in- 
creasing strand length. 

In polymer physics 0, El stresses are usually derived 
from chain conformations by dividing the chains into 
segments of length n. Segments are assumed to be- 
have as independent entropic springs with spring con- 
stant k B T 'jr 2 a (n). For a given monomer density p m , the 
segment density is p s {n) — p m /n [^J. Normal tensions 
derived from the virial tensor for this mesoscopic poly- 
mer model take a simple form involving only those ratios 
of mean-square internal distances which we have plotted 
in Fig. Eh: 



,(A) = k B T lim p s (n) 

n— >0 



•f(n,A||) 
r 2 a {n) 



(n, Aj 



r 2 a {n) 



(4) 

There are considerable subtleties in comparing the virial 
tensor calculated from the full microscopic interactions 
to the "Gaussian" normal tensions Fig. Eh shows 

that in the present case both quantities agree within the 
statistical error even though 

CTT .Gauss 

(A)/ct t (A) « 90%. 

In contradiction to the arguments put forward in Ref. jfj, 
we take the good agreement as quantitative evidence, 
that the tube represention of the network conformation 
in terms of linear paths through the network properly 
accounts for the relevant microscopic deformations. 

The final step of our analysis is the parameter-free com- 
parison in Fig. Eh of the measured stress-strain relations 
to the prediction of the ME theory 



ctt (A) 



(A?-l) fl (A,|)-(Al-l) fl (Ai) 



(5) 



iknT 



d 4 - 

T,A 



2 A a rf T ,s 



8 4(A Q ) d A TA + X 2 a d* 



T,B 



where we use the values of g?t,a and dr,B from our fits 
of the microscopic deformations. The measured and the 
inferred normal tensions agree within the statistical error. 
To demonstrate the importance of the proper treatment 
of the deformation dependence of the confinement, we 
have also calculated stress-strain curves from the WE 
theory dr,B = oo) and the HS- 

RP theory 0, EJ (d T {X a ) = *jKd T ^ d T , A = oo). 
The comparison in Figs. Eh an d b shows that attempts 
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FIG. 3: (color online) Mooney-Rivlin plot of normal tensions 
as a function of strain: normal tensions sampled using the 
virial tensor (large symbols as in Fig.^l, and "Gaussian" nor- 
mal tensions ^T, Gauss (small filled symbols) vs. theoretical 
predictions based on the analysis of the microscopic deforma- 
tions ( Fig. (a) WE crosslinker tube theory, (b) HS-RP 
entanglement tube theory, (c) ME theory. Symbols and lines 
have a 10% and 30% error, respectively. 

along these lines are restricted to estimates of the order 
of magnitude of the elastic response. 

To summarize, we have used computer simulations 
to determine strain-dependent localization, length scale 
dependent microscopic deformations and macroscopic 
stresses in model polymer networks. Fig. directly 
validates Edward's original tube concept: on a scaling 
level the same length scale characterizes the spatial lo- 
calization of chemical crosslinks and the crossover from 
local liquid-like to global solid- like behavior of the micro- 
scopic chain deformations. Closer inspection (Figs. QJEl 
c,d) reveals the different character of (and the crossover 
between) cross-link and entanglement dominated confine- 
ment anticipated by the double tube model. In par- 
ticular, Figs. El and show that this model can simul- 
taneously describe microscopic and macroscopic aspects 
of the response of our model networks to strain by ac- 
counting for the strain dependence of the (effective) tube 
diameter. While there clearly remain open questions 
concerning longitudinal fluctuations in the entanglement 
tube 0, |2j| , we would like to emphasize that our sim- 
ulations provide a substantially improved empirical ba- 
sis for addressing these and similar problems in the con- 
trolled development of statistical mechanical theories of 
rubber elasticity. Currently, we are analy sing the strain 
dependence of the primitive path mesh [29] in an at- 
tempt to systematically link the phenomenological tube 
model to the microscopic connectivity and topology of 
our model networks. Finally, we note that our simula- 
tions can also help to validate critical steps in the data 
analysis of (scattering) experiments addressing these is- 
sues. 
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